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Introduction. The wonderful book by Andrews, Askey, and Roy [5] is mainly 
devoted to special functions of hypergeometric type - to the plain and q-hypergeo- 
metric series and integrals. Shortly before its publication, examples of a third type 
of such functions, related to elliptic curves, began to appear. A systematic theory of 
elliptic hypergeometric functions was constructed in 2000-2004 over a short period 
of time. The present complement reviews briefly the status of this theory to the 
end of 2006. It repeats where possible the structure of the complemented book [2], 
and it is based mainly on the material |44j . 

The theory of quantum and classical completely integrable systems played a 
crucial role in the discovery of these new special functions. An elliptic extension 
of the terminating very well poised balanced g-hypergeometric series io¥>g with 
discrete values of parameters appeared for the first time in elliptic solutions of the 
Yang-Baxter equation [14] associated with the exactly solvable models of statistical 
mechanics [5]. The same terminating series with arbitrary parameters appeared 
in [47j as a particular solution of a pair of linear finite difference equations, the 
compatibility condition of which yields the most general known (1 + l)-dimensional 
nonlinear integrable chain analogous to the discrete time Toda chain. In [35j . 
the generalized gamma functions were investigated in detail, including one of the 
elliptic analogues of Euler's gamma function, which appeared implicitly already in 
Baxter's eight vertex model. The appearance of such mathematical objects was 
quite unexpected, since no handbook or textbook of special functions contained 
any hint of their existence. Some relatives of these functions were considered only 
in the old original papers by Barnes [6] and Jackson [19]. 

Generalized gamma functions. In the beginning of the XX-th century 
Barnes [6 suggested a multiple zeta function depending on m quasiperiods ujj G C: 

( m (s,u;uj)= ^2 ( I ry\ s ' & = n\u\ + . . . + n m Um, N = {0, 1, ...}, 

where the series converges for Re(s) > m and under the condition that if ujj /ujk G K, 
then ujj/uok > 0. Using an integral representation for £ m , Barnes also defined 
the multiple gamma function T m (u;ui) = exjp(dCm(s 1 u;uj)/ds)\ s= Q 1 which has the 
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infinite product representation 

1 k °° h 

— — =e Er=o7m.%- u tt ' c 1+ ") e Er= 1 (-i)^ ! (i) 

where jmk are some constants analogous to Euler's constant (in [BJ, the normal- 
ization 7 m o = was used) . The primed product means that the point n\ = . . . = 
n m = is excluded from it. The function T m (u;uj) satisfies m finite difference 
equations of the first order 

T m {u + Uj;u) =r^_ 1 (u;uj(j))T m (u;uj), j = l,...,m, (2) 

where ui(j) — (w\, . . . , wj-i, ^>j+i, ■ ■ ■ , w m ) and Tq(u;uj) := X/u. The function 
Ti(ojix;u)i) essentially coincides with Euler's gamma function T(x). The plain, q-, 
and elliptic hypergeometric functions are respectively connected to T m {u;u) for 
m = 1,2,3. 

We define three base variables p,q,r S C related to the pairwise incommensurate 
quasiperiods ^1,2, 3 as follows: 

27riii 27ri^ 2-n-i^ 

q = e "2, p = e "2 , r = e "1 , 

q = e "i , p = e "a , r = e "3 : 

where q,p,r denote the modular transformed bases. For < 1, the infinite 

products 

00 00 
(z;q)oc = - zq j ), {z;p,q) oc = ~[[(l-zp 3 q k ) 

j=0 j,k=0 

are well defined. It is easy to derive equalities [TH] 

= l-«, 7 r - = {z;P)oo, 7 — = (z;g) 00 . (3) 



{qz;q)ac " (qz;q,p)oo (pz;q,p) 

The odd Jacobi theta function (see formula (10.7.1) in [5]) can be written as 



i e 7rir(n+l/2) 2 e 7ri(2n+l)« 



6 1 (u\r) = -ij (-1) 
n= — 00 

= ^/ 8 e-™(p;p)^% 2 ™;p), ueC, 
where p = e 27 ™ T . The shortened theta function (see Theorem 10.4.1 in [2]) 

9{z-p) := (zMpz-^pU = — L- V(-l)V (fe - 1)/2 ^ fe (4) 

plays a crucial role in the following. It obeys the following properties: 

e(pz;p) = 0(z- 1 ;p) = -z- 1 9(z;p) (5) 
and 0(z;p) = for z = p k , k € Z. We denote 

0(oi, . . . ,a*.;p) := 6(ai;p) ■ ■ -6{ak\p), 0(at ±m ,p) := d(at;p)9(at~ 1 ;p). 
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Then the Riemann identity for products of four theta functions takes the form 

6{xw ± ,yz ± ;p) — 9(xz ± , yw^^p) = yw~ 1 9(xy ± , wz ± ;p) (6) 

(the ratio of the left- and right-hand sides is a bounded function of x G C* , and it 
does not depend on x due to the Liouville theorem, but for x — w the equality is 
evident). 

Euler's gamma function can be defined as a special meromorphic solution of 
the functional equation f(u + u>i) = uf(u). g-Gamma functions are connected to 
solutions of the equation f(u+Ui) = {l-e 27rlu ^ 2 )f(u) with q = e 2Al / W2 . For \q\ < 
1 , one of its solutions has the form T q (u) = 1 / ( e 27 ™'/"2 ; g) ^ defining the standard q- 
gamma function (it differs from function (10.3.3) in [5] by the substitution u = uj\x 
and some elementary multiplier). The modified g-gamma function ("the double 
sine" , "hyperbolic gamma function" , etc), which remains well defined even for \q\ = 
1, has the form 

/ f e ux dx\ 

7(U ; W ) = exp (- {1 _^ X){1 _^ X)T ) , (7) 

where the contour R + iO coincides with the real axis deformed to pass clockwise 
the point x = in an infinitesimal way. If Re(o>i), Re(w2) > 0, then the integral 
converges for < Re(u) < Re(wi + u>2)- Under appropriate restrictions upon u 
and 0*1,2, the integral can be computed as a convergent sum of the residues of 
poles in the upper half plane. When Im(wi/w2) > 0, this yields the expression 
j(u;uj) = (e 27rlu / UJl q;q) 00 /(e 2mu / UJ2 ;q) 00 , which can be extended analytically to the 
whole complex u-plane. This function, serving as a key building block of the ^-hy- 
per geometric functions for \q\ = 1, was not considered in [2] and [16]; for its detailed 
description see [501 1231 GUI [5D] and the literature cited therein. 

In an analogous way, elliptic gamma functions are connected to the equation 

f{u + u 1 )^e{e 2 ^^-p)f{u). (8) 

Using the factorization (0| and equalities (J3J) , it is not difficult to see that the ratio 

^{z-p,q) = — r (9) 

{z,P, Q)oc 

satisfies the equations 

F(qz;p,q) = 6{z;p)T(z;p,q), T(pz;p,q) = 9(z;q)T(z;p,q). 

Therefore the function f(u) = T(e 27Tm ^ 2 ;p,q) defines a solution of equation © 
valid for \q\, \p\ < 1, which is called the (standard) elliptic gamma function [35]. It 
can be defined uniquely as a meromorphic solution of three equations: ([H]) and 

f(u + w 2 ) = /(«), f(u + Ws ) = 9(e 2 ^/^ ; q)f(u) 

with the normalization /(X)L=i w ™/^) = •"■> smce there do not exist non-trivial 
triply periodic functions. The reflection formula for it has the form T(z;p, q) 
r(pq/z;p,q) = 1. For p = 0, we have r(z;0, q) = l/(z;q)oo. 
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The modified elliptic gamma function, which is well defined for |g| — 1, has the 
form [41] 

G(u- lu) = T(e 2 "^ ; p, q)T(re~ 2 ^ ■ q, r ). (10) 
It yields the unique solution of three equations: {8]| and 

/(« + wa) = e(e 2ffi «/ Wl ;r)/( U ), /(u + w 8 ) = e^^^/fa) 
with the normalization /(X)m=i = 1. Here 

„ , , U 2 U U LUl LU 2 1 

B 2l a(«;w) = + 7^ + 7T^ + o 

denotes the second Bernoulli polynomial appearing in the modular transformation 
law for the theta function 

#(V 2 ™^T;e- 2 ™Sr) = e ™ S2 ' 2 ^0(e 2 ™^;e 2 ™^) . (11) 

One can check |12j that the same three equations and normalization are satisfied 
by the function 

G(u;u) = e -™ p (">r(e- 2 ^;f,p), (12) 
where |f>|, \r\ < 1, and the polynomial of the third degree P(u) has the form 

^ 2 / 3ujiuj 2 lli 3 

■m—l / 

Functions (JTUJ) and (|T2"]) therefore coincide, and their equality defines one of the laws 
of the 5L(3; Z)-group of modular transformations for the elliptic gamma function 
[13] . From expression (|12|) , the function G(u;lu) is easily seen to remain meromor- 
phic when luxJlu^ > 0, i.e. \q\ = 1. The reflection formula for this function has 
the form G(a]Lu)G(b;u>) = 1, a + b = X^fc=i w k- In the regime \q\ < 1 and p, r — > 
(i.e., Im(aj3/wi), Im(w3/a'2) — > +00), expression (fT0)l obviously degenerates to the 
modified g-gamma function 7(1*; lo). Representation (|12p yields an alternative way 
of reduction to j(u; lu); a rigorous limiting connection of such a type was built for 
the first time in a different way by Ruijsenaars [35j . 

As shown by Barnes, the g-gamma function l/(z;q) 0o with z — e 2jrm /" 2 and 
q = e 27rlw i/"2 j \m{uii / LU2) > 0, equals in essence to the product F2(u;loi, 102)^2(11 — 
LU2 j lui 1 — LU2) • Analogously, the modified g-gamma function 7(11; to) coincides in es- 
sence with the ratio r 2 (wi +u>2 — u\ w)/r 2 (w,; lu). Since 9{z; q) = (z; q) o(<lz^ 1 ; q)oo, 
the r 2 (w; w)-function represents "a quarter" of the 9\{u/lu2\lui/lu2) Jacobi theta 
function (in the sense of the number of divisor points). Correspondingly, one can 
consider equation ||5J) as a composition of four equations for T3(u;lu) with different 
arguments and quasiperiods and represent the elliptic gamma functions as ratios 
of four Barnes gamma functions of the third order with some simple exponential 
multipliers [15] 141] . For some other important results for the generalized gamma 
functions, see [25| [30] . 
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The elliptic beta integral. It is convenient to use the compact notation 
r(ai, . . .,a k ;p,q) := r(en;p,g) • ■ • T(a k ;p, q), 

F{tz ± ',p, q) := T(tz;p, g)r(tz -1 ;p, q), T(z ±2 ;p, q) := T(z 2 ;p, q)T(z^ 2 ;p, q) 

for working with elliptic hypergeometric integrals. We start consideration from the 
elliptic beta integral discovered by the author in [38J . 

Theorem 1. We take eight complex parameters tj, j — 1, ... ,6, andp, q, satisfying 
the constraints \p\,\q\,\tj\ < 1 ondJ\j =i tj = PQ- Then the following equality is valid 

K Jr r(***; Ptg ) T = II rfet^ftg), (is) 

J 1 K iriij l<j<k<6 

where T denotes the positively oriented unit circle and k = (p;p) oa (q;q) 00 /4Tri. 

The first proof of this formula was based on the elliptic extension of Askey's 
method |3]. A rather short proof was given in [33], but still it does not fit the 
margins of this page and we skip it. The elliptic beta integral (p~3|) defines the 
most general known univariate exact integration formula generalizing Euler's beta 
integral. For p — > 0, one obtains the Rahman integral [5^ (see Theorem 10.8.2 in 
[!)), which goes to the well known Askey- Wilson g-beta integral [3] (see Theorem 
10.8.1 in [2]), if one of the parameters vanishes. 

We replace T by a contour C which separates sequences of poles converging to 
zero along the points z = tjq p m , k,m S N, from their reciprocals obtained by the 
change z — > 1/z, which diverge to infinity. This allows us to lift the constraints 
\tj\ < 1 without changing the right-hand side of (fl3|) . We substitute t$ — pq/A, 
A = nf=i*») and suppose that \t m \ < 1, m = 1, . . . ,4, \pt 5 \ < 1 < |t 5 |, < 
and that the arguments of t s , s — 1, . . . , 5, and p, q are linearly independent over 
Z. Then the following equality takes place [ID] : 

/" dz f dz 

k A E {z,t)— = K A E (z,t)— + c (t) V i/ n (t), (14) 

JC ^ ^ ^ |t 5 g«|>l,n>0 

where A E {z,t) = Y\ 5 m=1 T(t m z ± ;p,q)/r(z ±2 ,Az ± -,p,q) and 

„ (f) Ul=intmtt;p,g) , e{tjg 2n ;p) A 9(t m t 5 ) n 

C0( - j " r(t 5 - 2 ,A*± ;P , 9 ) ' Un[t) - e(t 2 ; P ) n o e {q t-%) n q ■ 

We have introduced here the new parameter t$ with the help of the relation 
rim=o — Q an d used the elliptic Pochhammer symbol 

o{t) n = n W;p) = r ff" :p 'f , . , := n 

(the indicated ratio of elliptic gamma functions defines 0(i) n for arbitrary n G C). 
The multiplier k is absent in the coefficient Co due to the relation lim 2 -*i(l — 



6 



V. P. SPIRIDONOV 



z)T(z;p,q) = l/(p;p)oo{q', q)oc and doubling of the number of residues because of 
the symmetry z — > z^ 1 . 

In the limit t^t^ — > q~ N , N £ N, the integral value in the left-hand side of (fT4")l . 
i.e. the right-hand side of (fT3f . and the multiplier co(t) in front of the sum of 
residues on the right-hand side diverge, whereas the integral over the unit circle 
T remains finite. After dividing all the terms by Co(t) and going to the limiting 
relation, we obtain the Frenkel-Turaev summation formula 

^ -' fit 1 9*5 9*5 9*5 \' v ' 

71=0 "V $1*2*3*5 ' *1 ' *2 ' *3 

which was established for the first time in [14] by a completely different method. 
For p — > and fixed parameters, equality (JT5J) reduces to the Jackson sum for a 
terminating 8</27-series (see formula (12.3.5) in [5]). 

General elliptic hypergeometric functions. Definitions of the general ellip- 
tic hypergeometric series and integrals were respectively given and investigated in 
detail in and [H]. So, formal series X^nez c « are called the elliptic hypergeo- 
metric series, if c n+1 = h(n)c n , where h(n) is some elliptic function of n € C. It 
is well known [5], that an arbitrary elliptic function h(u) of order s + 1 with the 
periods W2/0J1 and oj^/lj\ can be represented in the form 

w s tt 6(t k z;p) 

Hu) = y \\ —, r, (16) 

^Jl 9{Wkz;p) 

where z — q u . The equality h(u + uji/^x) — h{u) is evident, and the periodicity 
h(u+u>3/u)i) = h(u) brings in the balancing condition Ilfe=i tk — YVkt^i w k- Because 
of the factorization of h(n), in order to determine the coefficients c n it suffices to 
solve the equation a n+ i = 9(tq n ;p) a n , which leads to the elliptic Pochhammcr 
symbol a n = 9(i) n ciq. The explicit form of the bilateral elliptic hypergeometric 
series is now easily found to be 

( tl,...,t s +l \ TT ^(t^n n 

s+iG s+ i ; q,p; y := > M — — r- y n , 

\w u ...,w s+1 J ^^9W" 

where we have chosen the normalization Co — 1. By setting u> s +i — q, t s+ \ = to, 
we obtain the one sided series 

to,ti,...,t s \ s—^, 9(to,ti, . . . ,t s ) n 



\wi,...,w s ) ^9{q,w 1 ,...,w s ) n 

For fixed tj and wj, the function s +iE s degenerates in the limit p — > to the basic 
g-hypergeometric series s +i^ s with the condition Ilfc=o = lY[k=i w s- There 
are some problems with the convergence of infinite series (|17[) . and therefore we 
assume that they terminate due to the condition tk = q~ N p M for some k and 
N € N, M G Z. For consideration of some questions, the additive system of notation 
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for the elliptic hypergeometric series is more convenient (see, e.g., Ch. 11 in [16] or 
[44]), but we skip it here. 

Series (IT71) is called well poised, if t$q — t\Wi = . . . = t s w s . In this case the 
balancing condition takes the form t\ ■ ■ -t s = ±q^ s+1 ^/ 2 t^ ^ , and the functions 
h(u) and s +iE s become invariant under the changes tj — > ptj, j — 1, .... s — 1, 
and to — > p 2 to. For odd s and the balancing condition of the form t\- ■ -t s — 
-I^O+i)/ 2 ^ 1 )/ 2 j there appears the symmetry to — > pto and s +iE s becomes an 
elliptic function of all free parameters logiy, j — 0, . . . , s — 1, with equal periods 
(such functions were called in [391144] the totally elliptic functions). Under the four 
additional constraints i s _ 3 = t s -2 = -qVki, t s -i = qy/h/p, t s = -q\fph, 

connected with the doubling of the argument of theta functions, the series are called 
very well poised. In |40j , a special notation was introduced for the very well poised 
elliptic hypergeometric series: 



/ t ,h,..., i s _4, <7\/to, -qV^h Qy/to/p, -?VP*o 
\qt a /ti, qto/tg-i, v 7 ^, -\/to, \fpk,i~y'to/p 

~~m r" 11 at + ( qy > =: s+iV s (t ;t 1 ,...,t s _ 4 ;q,p;y), 

where the balancing condition has the form rifc=i ^fe = 5 ^ 2 q( s ~ 7 ^ 2 , and for 

odd s we assume the positive sign choice for preserving the symmetry to — » pto. For 
the y argument value y = 1, it is omitted in the series notation. Summation formula 
(|15p thus gives a closed form expression for the terminating ioVg(io; t±, . . . , t§; q,p)- 
series. 

Contour integrals J c A(u)du are called the elliptic hypergeometric integrals, if 
the function A(u) satisfies the system of three equations 

A(u + ij k ) = h k (u)A(u), k= 1,2,3, (19) 

where LUi.2,3 G C are some pairwise incommensurate parameters and hk(u) - some 
elliptic functions with periods cut, ^k+i (we set ujk+3 — Uk)- One can weaken 
requirement (j 1 Q j) by leaving only one equation, but then there appears a functional 
freedom in the choice of A(it), which should be fixed in some other way. 

Omitting the details of consideration [5TJ 03] , we present the general form of 
permissible A(u). We suppose that this function satisfies equations (fl9j) for fc = 1,2, 
where 

hl( " )=Itt !]9( Wj e««/«»;p)' ^W=y 2 ll^. e _ 27rWwi;r) , 
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with \pl\r\ < 1 and ]T*=i *j = Hj=i w j> IIj=i *j = 11^=1 If we take M < *> 
then the most general meromorphic function A(u) has the form 

s T-i/-, 27Ti— \ ^ T-i/r -27ri— - \ m ni 2ni— \ 

aw = n rfee 2 : ;p ' g) n fce 2 : ;g,r) tt (afc£ 2 ^\ m 

j=i r ( w je -2 ; p, g ) j=1 r(m,e <"i ; q, r) fe=1 6<(& fe e "2 ; g ) 

where the parameters d and m are arbitrary, and dk,bk,c are connected with y\ 
and 2/2 by the relations t/2 = e c " 2 and j/i = e caJl ri/cLi ^fc a fe • ^ appears that the 
function h^{u) cannot be arbitrary - it is determined from expression f|20[) . 

For \q\ = 1, it is necessary in (j2"0|) to choose I = s and to fix parameters in such 
a way that the T-functions are combined to the modified elliptic gamma function 
G(u\oS) (it is precisely in this way that this function was built in [41]): 

A(u = G(u + g s ; U ) eCU+d) 

where the parameters gj , Vj are connected with tj , Wj by the relations tj — e 27ri9j /" 2 , 
wj = e^ivj/^^ an( j y i 2 _ e cwi,2 rpjjg integrals J c A(u)du with kernels of the in- 
dicated form define elliptic analogues of the Meijer function. A more general theta 
hypergeometric analogue of the Meijer function was constructed in [41j . but we skip 
it here. 

We limit consideration to the choice t = m = in (12TJ)) . The correspond- 
ing integrals are called well poised, if iiioi = . . . = t s w s = pq. The additional 
condition of very well poisedness fixes eight parameters t s — 7, . . . , t s — {±(p<7) 1//2 , 
±q 1 / 2 p : ±p 1 / 2 q, ±pq} and doubles the argument of the elliptic gamma function: 
r(tj2;p, 9) = l/r(z _2 ;p, 9). The most interesting are the very well poised 
elliptic hypergeometric integrals with even number of parameters 

K"> • ■ • , h m+e ) = k [ u % T ^\ pt q) * , n +6 ^w ro+i , (22) 

Jt T(z ±2 ;p,g) 2 

with \tj\ < 1 and the "correct" choice of the sign in the balancing condition. They 
represent integral analogues of the s+ iVs-series with odd s, "correct" balancing 
condition and the argument y — 1, in the sense that these series appear as residue 
sums of particular pole sequences of the kernel of I^ m \ We note that 1^ coincides 
with the elliptic beta integral. 

Properties of the elliptic functions explain the origins of the old hypergeometric 
notions of balancing, well poisedness, and very well poisedness. However, strictly 
speaking these notions are consistently defined only at the elliptic level, because 
there are limits to such q-hypergeometric identities in which they are not conserved 
any more [551 HO]! The fact of unique determination of the balancing condition 
for series (TT5|) with odd s and integrals (|2"2")l (precisely these objects emerge in the 
main part of interesting applications) illustrates a deep internal tie between the 
"elliptic" and "hypergeometric" classes of special functions. Multivariable elliptic 
hypergeometric series and integrals are defined analogously to the univariate case 
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V(ti,...,ts;p,q) = k I ' — ( 23 ) 



- it is necessary to use systems of finite difference equations for kernels with the 
coefficients which are elliptic functions of all summation or integration variables 
[391 141] , which is a natural generalization of the approach of Pochhammer and 
Horn to the functions of hypergeometric type [2 [IT] . 

An elliptic analogue of the Gauss hypergeometric function. We take 
eight parameters t\, . . . , tg £ C and two base variables p,q £ C satisfying the con- 
straints \p\, \q\ < 1 and Yij=i tj ~ P 2( l 2 (^ ne balancing condition). For \tj\ < 1, j = 
f , . . . , 8, an elliptic analogue of the Gauss hypergeometric function 2-^1(0, b; c; x) is 
defined by the integral representation [44] 

Y\]=i^^z ± ;p 1 q) dz 
T(z ±2 ;p,q) z 

i.e. by the choice m = 1 in expression (|22[) . For other admissible values of pa- 
rameters, the ^-function is defined by the analytical continuation of expression 
(|23p . From this continuation one can see that the V- function is meromorphic in all 
parameters for arbitrary tj G C* — it is sufficient for this to compute residues of 
the integrand and to define analytically continued function as a sum of the integral 
over some fixed contour and residues of the poles crossing this contour. For a more 
detailed analysis of the meromorphic character of the elliptic hypergeometric inte- 
grals and an interesting role of the "correct" sign choice in the balancing condition 
played in that, see [2"7| . 

The first nontrivial property of function (|23D consists in its reduction to the ellip- 
tic beta integral under the condition tjt^ = pq for an arbitrary pair of parameters 
tj and tfc (expression (|f 3p appears from t^tg = pq). The ^-function is evidently 
symmetric in p and q. ft is invariant also under the S^-group of permutations of 
parameters tj isomorphic to the A7 Weyl group. We consider the double integral 

f ITLi F(ajZ ± ,bjW ± ;p, q) T(cz ± w ± ;p, q) dz dw 
J T 2 T(z ±2 ,w ±2 ;p,q) z w ' 

where Oj, bj, c £ C, \a,j\, \bj\, |c| < f , and c 2 Ytj=i a j — ^ bj = pq- Using for- 
mula (|13|) for integrations over z ot w (the permutation of the order of integrations 
is permitted), we obtain the following transformation formula: 

V(t)= Y[ T(t J t k ,t J+4 t k+4 ;p,q)V{s), (24) 

l<j<fc<4 

where we used the compact notation V(t) = V{t\, . . . ,ts;p, q), 

Sj=P~ 1 tj, j = 1,2,3,4 _ jtitjtiU _ j pq 



s j = Ptj, j = 5,6,7,8 ' 9 V Pq V hht 7 t% 

and |tj|,|s 3 '| < 1. This fundamentally important relation was obtained by the 
author in [41] , where the function V(t) appeared for the first time, ft represents an 
elliptic analogue (moreover, its integral generalization) of the Bailey transformation 
for four non terminating io(/3g-series |16j . 
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We repeat transformation (|24[) once more with the parameters 53,4,5,6, playing 
the role of £1,2,3,4, and permute parameters £3, £4 with t5,tg in the result. This 
yields the relation 

4 

vit) = [] r(tjt k+4 ; P , q) v(Ti/ti, r*/t 4 , u*/t 6 , u?/t 8 ), (25) 

3,k=l 

where T = t^hU, U = t 5 t 6 t 7 t$ and iTp/ 2 < \tj\ < 1, |C/| 1/2 < \t j+4 \ < 1, j = 
1,2,3,4. We equate now the right-hand sides of relations (pM)) and (f2"5j) , express 
parameters tj in terms of Sj and obtain the third relation 

V(s)= Y[ T ( s 3 s k;P,q)V(^/pq/s 1 ,...,^/pq/s 8 ), (26) 

l<j<fe<8 

where ImI 1 / 2 < \sj\ < 1 for all j. 

We consider Euclidean space M 8 with the scalar product (x, y) and an orthonor- 
mal basis e\ £ K 8 , (ei,ej) = Sij. The root system A 7 consists of the vectors 
v = {ei — ej, i 7^ j}. Its Weyl group is composed from the reflections x — ► S v (x) — 
x — 2v(v, x)/(v, v ) acting in the hyperplane orthogonal to the vector X)i=i e « (i- e -; 
the coordinates of the vectors x — Y2i=i x i e i satisfy the constraint Y^=i x i = 0)> 
and it coincides with the permutation group S%. 

We connect parameters of the ^(i)-function with the coordinates Xj by the 
relations tj = e 2vlXj (pq) 1 ^ 4 . Then the balancing condition assumes that Y^l=i x * = 
0. The first coordinate transformation (|24p is now easily seen to correspond to 
the reflection S v (x) for the vector v = QZ*_ 5 e » — Si=i e i)/2 having the canonical 
length (v,v) — 2. This reflection extends the group A 7 to the exceptional Weyl 
group £7. Relations (|25|) and ([26]) were proved in a different way by Rains in [27] . 
where it was indicated that these transformations belong to the group Ej. 

We denote by V(qtj,q~ 1 tk) elliptic hypergeometric functions contiguous to V(t) 
in the sense that tj and t k are respectively replaced by qtj and q~ 1 t k . The following 
contiguous relation for the V^-functions is valid 

t 7 6 (t s 4/q;p) V(qt 6) q-H s ) - (t 6 t 7 ) - t 7 (t 6 t^;p) V(t), (27) 

where (tg i 7 ) denotes the permutation of parameters in the preceding expres- 
sion (such a relation was used already in [38]). Indeed, for y — te,w — t 7 , and 
x = q~ 1 t% the Riemann identity §E§ is equivalent to the q-difference equation for V- 
function's integrand A(z,t) = Yl k=1 T(tkZ ± ;p, q)/T(z ±2 ;p, q) coinciding with P7|) 
after replacement of V"-functions by A(z,i) with appropriate parameters. Integra- 
tion of this equation over the contour T yields (|27|) . We now substitute symmetry 
transformation (|26| in (|27jl and obtain the second contiguous relation 

5 5 

ted {T s ;p ) II »(— ;p)^(9 _1 *8»9*8) - (*e - *t) = n «(*8**;p))V(t). 

^ ^ fc— 1 ^ ^ A.:— 1 
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An appropriate combination of these two equalities yields the equation 

A{t){u{qt 6 ,q-H 7 ) - U(tj) + fa «- t 7 ) + U(t) = 0, (28) 
where we have denoted U(t) — V(t)/T(tQtf,t7t%;p,q) and 

U e(t 6 /t7,t 7 /qt 6 ,t 6 t 7 /q;p) W 0(t 8 t fc ;p) 

Substituting t,- = e 2 ^ l9j ^ 2 , one can check that the potential _4(i) is a modular 
invariant elliptic function of the variables gi, ■ ■ ■ ,37, i.e. it does not change after 
the replacements gj — > + 0^2,3 and (w2,W3) (— 1^3, ^2)- 

We introduce now the new notation for parameters: t§ — cx, t 7 — c/x, and 

£fe = — , k = 1, . ..,5, £ 8 = —, £7 = — ■ 
ct k t 8 q 

Setting c = ^Js^e 8 jp 2 \ we find that the balancing condition in terms of e k takes 

the form Ilfe=i £ k — P 2( l 2 ■ After the replacement of U(t) in (|28[) by some unknown 
function f(x), we obtain a (/-difference equation of the second order which is called 
the elliptic hyper geometric equation: 

A(x) (f(qx) - f(x)) + Aix- 1 ) (fiq-'x) - f(x)) + uf(x) = 0, (30) 



k=l 



We have already one functional solution of this equation 

= V{q/c£i, . . . ,q/ce 5 ,cx,c/x,c/e 8 ;p, q) ^ 
T(c 2 x ± /e 8 ,x ± e 8 ;p,q) 

where, for the choice C = T of the integration contour in the definition of the V- 
function, it is necessary to impose the constraints (in the previous parametrization) 
y/\pq\ < \tj\ < 1, j = 1,...,5, and y/\pq\ < l^^el, \q ±X t 7 \, k^sl < 1, which can 
be relaxed by analytical continuation. Other independent solutions can be obtained 
by the multiplication of one of the parameters £1, . . . , £5, and x by some powers of 
p or by permutations of £1, . . . , £5 with eg. 

We denote e k = e 27rWfc /" 2 , x = e 2 ™/" 2 , and f x {x) =: Fi(u; a; ojx, uj 2 , a; 3 ). After 
these changes one can check that equation (|30l) is invariant with respect to the 
modular transformation (^2,^3) — > (—^3,^*2). One of the linear independent solu- 
tions of (|30"f therefore has the form F 2 (u; a; u>i, u> 2 , 0J3) := Fi(u; a; u>i, — LJ3, 102)- The 
same solution would be obtained if we repeat the derivation of equation (|30p and its 
solution (|32[) after replacing the T-functions by the modified elliptic gamma func- 
tion G{u;(jj). The ^-function is therefore well defined even for \q\ = 1. Different 
limiting transitions from the V^-function and other elliptic hypergeometric integrals 
to g-hypergeometric integrals of the Mellin-Barnes or Euler type are described in 
[44), [45] and much more systematically in [8| [30] . 
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Biorthogonal functions of the hypergeometric type. In analogy with the 
residue calculus for the elliptic beta integral (fT"4|) . one can consider the sum of 
residues for a particular geometric progression of poles of the ^-function kernel for 
one of the parameters. This leads to the very well poised 12 Vii-elliptic hypergeomet- 
ric series the termination condition of which is guaranteed by a special discretization 
of the chosen parameter. In this way one can derive contiguous relations for the 
terminating i^V^i-series [47l [48] out of the contiguous relations for the ^-function, 
which we omit due to the lack of space. For instance, this yields the following 
particular solution of the elliptic hypergeometric equation (|30[) : 

u r \ t/ ( £ s 1 1 1 IP IP £ 6 \ /oca 

Rn{x;q,p) = 12V11 — ; , , , , ,£52:, — ;q,p , (33) 

\ES £1^8 £2^8 £3^8 £4^8 £5£s x / 

valid under the condition pq/e^es = q~ n , n € N (we remind that Ilfc=i £ k — 
p 2 q 2 )- We describe below properties of the i? n - function found in [31] in the notation 
passing to our after the replacements to, 1,2 — > £1,2,3, £3 — + £6, £4 ~* £s, A 1 - * £4£s/.P<7, 
and Afi/qti — > pq/e^es- 

Equation (130|) is symmetric in £1,.. ■,£&■ Since series (|18|l are elliptic in all pa- 
rameters, function (|33|> is symmetric in £1, . . . , £5 and each of them can be used for 
the termination of series. A permutation of £1, £2, £3, £5 with £6 yields R n {z; q,p) up 
to some multiplier independent on x due to an elliptic analogue of the Bailey trans- 
formation for terminating i2Vn-series |14) . which can be obtained by degeneration 
from equality (pM]) . 

The same contiguous relations for the 12 Vn -series yield the following three term 
recurrence relation for R n (x; q,p) in the index n: 

(z(x) ~ a n+ i)p{Aq n ^ 1 /e 8 ) (R n+ i(x; q,p) - R n {x; q,p)) + (z(x) - /3 n -i) (34) 
xp(q~ n ) (R n -i(x;q,p) - R n (x;q,p)) + 5(z(x) - z(s 6 ))R n (x; q,p) = 0, 

z{x) = d a f^% a n = z(q n /e 8 ), f3 n = z{Aq n -% 



P(t) = 



2 + 

£6 g£6 qt qt qt q trj 
Sst ' £8* ' £1£2 ' £2£3 ' £1£3 ' A 



:p 



( 



qt 2 e 8 q 2 t 2 eg . 
A ' A ■ 



g 2 £e q q q ± 
—7-, — , — , — ,£e?y ;p 

A £i£8 £2£8 £3£8 

where A = £i£2£3£6£s, and £ and r\ are arbitrary gauge parameters, £ 7^ r] ±1 p k , k S 
Z. The initial conditions = and i? = 1 guarantee that all the dependence 
on the variable x enters only through z(x), and that R n {x) is a rational function 
of z{x) with poles at the points a%, . . . , a n . 

The elliptic hypergeometric equation for the i?„-function can be rewritten in the 
form of a generalized eigenvalue problem T>\R n = \ n T>2R n for some g-difference 
operators of the second order T>\^ and discrete spectrum X n [41]. We denote by <f>\ 
solutions of an abstract spectral problem T>\<j)\ — XT>2(f>x, and by solutions of the 
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equation T>fip\ — XD^ipx, where T>i 2 are the operators conjugated with respect 
to some inner product (fp\(j)}, i.e. (Pi^^) = (^1^1,20). Then = (VvK^i — 
XT>2)(j)\) = {n — X){P2'4'n\<j>\) 1 i.e. the function X^Vv is orthogonal to for 
/i 7^ A. This simple observation was put by Zhcdanov into the basis of the theory of 
biorthogonal rational functions generalizing orthogonal polynomials [53] (see also 
[J71 HE]). Analogues of the functions l^VVi for R n (z; q,p) have the form 

t I \ f Ae 6 A A A e 6 qp qp \ , . 

T n (x]q,p) = 12V11 ; — , — , — ,E 6 x, — , , ;q : p , (35) 

V q £1 £2 £3 x £ 4 £ 8 £5£8 / 

which are rational functions of z{x) with poles at the points /3i, . . . , [3 n . 

We denote R nm {x) = R n (x;q,p)R m (x;p,q) and T nm (x) = T n (x; q,p)T m (x;p, q), 
where all the i2Vn-series terminate simultaneously because of the modified termi- 
nation condition £ 4 £ 8 = p m + 1 q n + 1 ^ n,m G N. The functions R nm solve now not 
one but two generalized eigenvalue problems which differ from each other by the 
permutation of p and q. 

Theorem 2. The following two-index biorthogonality relation is valid: 
r ( M? 1 ^lhes^j^jP^)dx 

T n i(x)R mk (x) , T r = h n i d mn ku (36) 

n ,ki T(x ±z ,Ax ± ;p,q) x 

where S = {1,2,3,6,8}, C rnn ja denotes the contour separating sequences of points 
£j p a q b (j = 1,2,3,6), e 8 p a - k q b - m ,p a+1 - l q b+1 - n / A, a, b S N, from their x -> x" 1 
reciprocals, and the normalization constants have the form 

, rij<fc,j,fcg5r(£jg fc ;p, g ) 

«nz = ' ' , . _! : — K(q,P)-hi(p,q), 



h n {q,p) 



8(A/qe 8 ; p)0{q, g£e/£s, £i£2, £i£3, £2£3, ^£e)n g n 
6»(A(7 2n /q£8;p)^(l/£6£8, £l£6, £2£6, £3£6, A/qe 6 , A/qe$) n ' 



This theorem was proved in [41] by direct computation of the integral in the left- 
hand side with the help of formula (fl"3]) . Appearance of the two-index orthogonality 
relations for functions of one variable is a new phenomenon in the theory of special 
functions. It should be remarked that only for k = I = there exists the limit 
p — > and functions R n (x; q,0), T n (x; q,0), appearing from that, coincide with 
Rahman's family of continuous ioyg-biorthogonal rational functions [26 . A special 
limit Im(aj3) — > 00 in the modular transformed R nm and T nm leads to the two-index 
biorthogonal functions which are expressed as products of two modular conjugated 
io</>9-series [44] . A special restriction for one of the parameters in R n (x]q,p) and 
T n (x;q,p) leads to the biorthognal rational functions of a discrete argument [47] 
generalizing Wilson's functions [S3]. All these functions are natural generalizations 
of the Askey- Wilson polynomials |4J. Note that R nm (x) and T nm (x) are mero- 
morphic functions of x € C* with essential singularities at x — 0, 00 and only for 
fc = ( = 0orii = m = 0do they become rational functions of some argument 
depending on x. It is expected that there are more general systems of biorthogonal 
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functions based on the T^-function and connected to solutions of the generalized 
eigenvalue problem with continuous spectrum or recurrence relation (]34|) with the 
shifted index n. 

Elliptic beta integrals on root systems. We introduce an analogue of the 
constant k for the C n (or BC n ) root system K n = (p;p)^ (q;q)2 /(27ri) n 2 n n\ and 
describe a C„-elliptic beta integral representing a multiparameter generalization of 
integral (fT3| . 

Theorem 3. We take n variables Z\,...,Z n € T and complex parameters t%,..., 
t2n+4 andp,q satisfying the constraints \p\, \q\, \tj\ < 1 and Hj=t 4 tj = PQ- Then 

tt 1 yr Tlm=ir(t m zf;p,q) dzj dz n 

T(z?z±;p,q) /j T(zf;p,q) 



11 T(t m t s ;p,q). (37) 



l<m<s<2n+i 



Formula ([37]) was suggested and partially confirmed in [TT], and it was proved 
by different methods in [27] [43] [44] . It reduces to one of Gustafson's integration 
formulas [18] in a special p — > limit. 

Theorem 4. We take complex parameters t,t m (m = 1, . . . ,6),p and q restricted 
by the conditions \p\,\q\, \t\, \t m \ < 1, and t 2n ~ 2 Y[ m= itm = pq- Then 

tt T (tzfzf;p,q) " Y\l l=1 T{t m zf\p,g) dzi dz n 
1</<L nzfzt,P,9) M F(zf ;p, ff ) * " " * *, 



r(t J ;p,g) 

, ' , ' : ''-'/' 1 ,„ ,J, 



H r(ti-H m t s ;p,q) . (38) 



In order to prove ([55]) , we consider the following (2n — l)-tuple integral 



n i A n^pTftr ^ ; p> q) 



^ 



1 



n nt^zfw^q) n 



i<i< fe <n-i r K ±u, fc;P'9) 

l<fc<n-l 

TT r («^* n ~ 3/2 n"=i*.;P.9) dwi dwn-xdzx dz n 
x s > (39) 

with g, £ and i r , r = 0, . . . , 5, lying inside the unit circle such that t n ~ [\r=o tr — 
pq. We denote the integral in the left-hand side of ([38]) as I n (t, t\, . . . , t 5 ;p, q). 
Integration over the w-variables with the help of formula (|37| brings expression 
([39]) to the form r 7l (i)r -1 (i") I n (t, t%, . . . , ts;p, q), where we have introduced the 
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parameter t& via the relation t 2n ~ 2 Ylj=i tj = PQ- Because the integrand is bounded 
on the integration contour, we can change the order of integrations. As a result, 
integration over the ^-variables with the help of formula ([57]) brings expression (j3T))) 
to the form T n -\t) n o <r< s < 5 r (^* s ) i»-i(t, t 1/2 *i, ■ ■ ■ , t^H^p, q), i.e. we obtain 
the following recurrence relation in the dimensionality of the integral of interest n: 

I n {t,t u ...,t 5 ;p,q) = iP 'f IT r (*^ S ;P, 9) /n-l (*, t 1/2 h , . . . , i 1 /^ I P, q). 

0<r<s<5 



Iterating it with known initial condition (|13[) for n = 1, we obtain f|38[) . 

Integral (|38[) was suggested by van Diejen and the author in [10] . The given proof 
is taken from [11] , it models Anderson's proof of the Selberg integral described in 
[2] (see Theorem 8.1.1 and Sect. 8.4). It represents also a direct generalization of 
Gustafson's method [TB] of derivation of the multiple g-beta integral obtained from 
(|38p after expressing t§ via other parameters, removing the multipliers pq with the 
help of the reflection formula for T(z;p, q), and taking the limit p — > 0. A num- 
ber of further limiting relations in parameters leads to the Selberg integral - the 
fundamentally important integral because of many applications in mathematical 
physics. Therefore formula (|3"5]) represents an elliptic analogue of the Selberg in- 
tegral (an analogous extension of Aomoto's integral described in Theorem 8.1.2 of 
[2] is derived in [27] ) . It can be interpreted also as an elliptic extension of the BC n 
Macdonald-Morris constant term identities. 

There are several other elliptic beta integrals on root systems. In particular, the 
author has suggested three different integrals for the A n root system (two of them 
look different for even and odd values of n) in [?T] . In [13] , Warnaar and the author 
have found one more A„-integral which appeared to be new even after degeneration 
to the q- and plain hypergeometric levels. 

In analogy with the one dimensional case [41] . it is natural to expect that the 
multiple elliptic beta integrals define measures in the biorthogonality relations for 
some functions of many variables generalizing relations ([3T>|) . In [27J [2H], Rains 
has constructed the first system of such functions on the basis of integral ([38]) . 
These functions generalize also the Macdonald and Koornwinder orthogonal poly- 
nomials, as well as the interpolating polynomials of Okounkov. In this sense, the 
results of [271 EH] represent to the present moment the top level achievements of 
the theory of elliptic hypergeometric functions of many variables. The author espe- 
cially likes the following £?C n -generalization of transformation ([2"4"| proven in [27J : 
I(t 1 ,...,t 8 ;t;q,p) = I(si, . . . , s 8 ; t; q,p), where 

I(t 1 ,...,t 8 ;t;q,p) f yr T{tzfz±;p,q) " ULi^zf ; p, q) dz, 

ni<,< fe <8 r (^;^9,t) Kn k^4} k<n n4*frp, q) M r(z± 2 ;p, g ) zj : 



ptj, j = 5,6,7,8 1 ^ V Pqt X ~ n V*5*6*7*8' 



= P~H j , i = 1,2,3,4 . _ ft^_Jpqt^ 1 \ tmAsjl<1: 
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and T(z:p,q,t) = H°° k:l=0 (l ~ zPp k q l ){\ - z- 1 V +1 p k+1 q l+1 ) is the elliptic gamma 
function of the higher level connected to the Barnes gamma function Ti(u;u>). 

Integral analogues of the Bailey chains. The Bailey chains, discovered by 
Andrews, serve as a powerful tool for building constructive identities for hypergeo- 
metric series (see Ch. 12 in [2]). They describe mappings of given sequences of 
numbers to other sequences with the help of matrices admitting explicit inversions. 
So, the most general Bailey chain for the univariate g-hypergeometric series sug- 
gested in [1] is connected with the matrix attached to the sV?7 Jackson sum [7]. An 
elliptic generalization of this chain for the s+ iVs-series was built in [40], but we do 
not consider here neither it nor its complement [52] • Instead of that, we present a 
generalization of the formalism of Bailey chains to the level of integrals discovered 
in 02]. 

We say that two functions a(z, t) and /3(z, t) form an elliptic integral Bailey pair 
with respect to the parameter t, if they are connected by the relation 

/3(w,t)=K / T(tw ± Z ± ;p,q)a(z,t) — . (40) 
Jt z 

Theorem 5. For a given elliptic integral Bailey pair a(z, t) and (3(z, t) with respect 
to the parameter t, the functions 

T(tuw ± ;p, q) 
T(ts 2 uw ± ;p, q) 



a {w,st) = j -a{w,t), 



T(t 2 s 2 ,t 2 suw ± ;p,q) f T(sw ± x ± ,ux ± ;p, q) dx 



r(s 2 , t 2 , suw^'jP, q) Jf T(x ±2 ,t 2 s 2 ux^;p, q) ' x 

where w € T, form a new Bailey pair with respect to the parameter st, and the 
functions 

. T(s 2 t 2 ,uw ± ;p,q) [ T{i?sux ± ,sw x^Pjq) . .dx 

a (w,t) = k — , t r / r r a(x,st) — , 

1 {s z , t z 1 , t^s^uw* ;p, q) Jj 1 [sux^ip, q) x 

T(tuw ± ;p,q) 

P [W, t) = — j r(3(w, St) 

1 (ts^uw^'jP, q) 
form a new Bailey pair with respect to the parameter t. 

The proof is sufficiently simple. In the first case, it is necessary to substitute 
the key relation for f3(x,t) in the definition of j3'(w,st), to change the order of 
integrations and to take off one of the integrations with the help of the elliptic 
beta integral (we omit restrictions for the parameters necessary for validity of this 
procedure). The second relation is proved in a similar way. These chain substitution 
rules introduce two new parameters u and s at each step of their application. In 
fact, they are related to each other by the inversion of the integral operator entering 
the definition of integral Bailey pairs [13] . 

This theorem is used in a way analogous to the scries case: it is necessary to take 
initial a(z,t) and (3(z,t), found, say, from formula (|13|) . and to generate new pairs 
with the help of the described chain substitution rules. Equality ([40l for these pairs 
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leads to a binary tree of identities for elliptic hypergeometric integrals of different 
multiplicities. As an illustration, we would like to give one nontrivial relation. With 
the help of formula (|13[) . one can easily verify validity of the following recurrence 
relation 

, n 2m+5<k<l<2m+S 

T(t k ti;p,q) 

1 K '{ti, ... ,t 2m +s) - ; (41) 

nl=Zm+5 T (Pm t k w ± ;p,q) (m) dw 

I y >{ti,...,t2m+A,Pm,W,p m W ) , 



T(w ±2 ;p, q) 

where p 2 n — rifc=2m+5 an ^ the integral /( m ) was defined in (22). By an 

appropriate change of notation for parameters, we obtain a concrete realization of 
the Bailey pairs: a oc I 1 *™ 1 * 1 and (3 cx /( m+1 ). For m = 0, substitution of the explicit 
expression for dHJ) in the right-hand side of flU]) yields identity tpij). Other 
interesting consequences of recursion (141) (an elliptic analogue of formula (2.2.2) 
in [2]) are considered in [44]. Various generalizations of the integral transformation 
(|40|) to root systems and their inversions are described in [46] . 

Conclusion. We would like to finish by listing some other achievements of the 
theory of elliptic hypergeometric functions. Multiple elliptic hypergeometric series 
were considered for the first time by Warnaar [51) . We described mostly properties 
of the elliptic hypergeometric integrals, since many results for the elliptic hyper- 
geometric series represent their particular subcases and they can be derived via 
the residue calculus. A combinatorial derivation of the Frenkel-Turaev summation 
formula is given in [37j . Various generalizations of this sum to the root systems were 
found in [TUJ [3TJ STJ [!)T] , and multivariable analogues of the Bailey transformation 
were described in [211 H3 ED IE] ■ Expansions in partial fractions of the products of 
theta functions and the identities connected to them were discussed by Rosengren 
[3"T] (see also [HJEZ])- The terminating continued fraction generated by three term 
recurrence relation (|34|) and the Racah type termination condition was computed 
in 48 . The raising and lowering operators connected with rational functions were 
discussed in (27] [49] . Connection to the Sklyanin algebra is considered in [24] [28] 
I32| 133] . In particular, in [33] Rosengren proved an old Sklyanin conjecture on the 
reproducing kernel. A systematic treatment of the elliptic determinant formulas 
connected to the root systems is given in the work of Rosengren and Schlosser [34] . 
In [55], it was shown that the i2Vn-series appears as a particular solution of the 
elliptic Painleve equation discovered by Sakai [36]. An analogous role is played 
by the general solution of the elliptic hypergeometric equation [44] [45] and some 
multiple elliptic hypergeometric integrals [29] . The page limits of the present review 
did not allow the author to cite a number of other interesting results, an essentially 
more complete review of the literature is given in [44] . 

Derivation of the i2Vn-function from a similarity reduction of some discrete 
integrable system, given in [47], reflects the essence of a heuristic approach to all 
special functions of one variable. This approach is described in detail in [44] using 
a number of examples, including some other new special functions. 
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The main part of the plain hypergeometric constructions admits thus natural 
elliptic generalization, although such parallels emerge for a rather large number 
of free parameters and structural restrictions. Nevertheless there remain many 
open problems, in particular, analysis of the conditions of convergence of infinite 
elliptic hypergeometric series, investigation of specific properties of the points of 
finite order on the elliptic curve, computation of the nonterminating elliptic hy- 
pergeometric continued fraction, detailed analysis of the non-self-dual biorthogonal 
functions of [47) and construction of their multivariable analogues, and so on. 

I am indebted to Yu. A. Neretin for the suggestion to write this complement and 
to G. E. Andrews, R. Askey, and R. Roy for an enthusiastic support of this idea. I 
am grateful to the Kyoto University (Research Institute for Mathematical Sciences 
and Graduate School of Informatics) for a hospitality and a support of my stay, 
during which the material of this review was presented in a lecture series. This 
work is partially supported by the Russian foundation for basic research (RFBR), 
grant no. 05-01-01086. 
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